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ABSTRACT 



An analytical model of a jet injected normally from a 
flat plate into a uniform crossing flow was modified to 
provide an improved method of predicting the interference 
effects arising from the complex flow fields induced by ship 
bow-thrusters. This model is an extension of previous work 
based upon a description of the jet as a series of distrib- 
uted vortices. The analysis takes into account the position 
of the effective source of the jet and the blockage due to 
the presence of the jet in the crossflow. Improvement in 
the predicted pressure distribution on the plate was obtain- 
ed by adopting a new formula for determining the effective 
jet source which more closely approximates the experimental 
results. Further improvement resulted from the incorpora- 
tion of a symmetric foil shape to simulate the presence of 
the jet and its highly turbulent wake. The accuracy of the 
model was evaluated for representative j et-to-crossf low 
velocity ratios. Good agreement with experimental results 
was achieved for large portions of the interaction field 
with particular improvement in the near-field and down- 
stream regions. 
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I. 



INTRODUCTION 



Since man first set sail, ships and boats have been 
steered by stern-mounted rudders or rudder-like devices. 

The bow has been left relatively uncontrolled. As the size 
of ships increased, they became more difficult to maneuver 
in close quarters, as when entering port. Tug boats were 
often required to provide lateral thrust at both ends of the 
ship. This was a satisfactory solution to the problem for 
many years. It did not provide the degree of independence 
desired by ship operators, however, and other methods were 
considered. Although experimental lateral thrust units were 
used in the mid-nineteenth century, the first practical unit 
was not installed until 1955. This unit, a through- tunnel 
type bow thruster installed in the "Princess of Vancouver" 
proved highly successful [1] . Since then, bow thrusters 
have steadily gained in popularity. The advantages of these 
units include providing a greater degree of maneuverability, 
increased safety, increased ability to operate in rough 
weather at oil rigs and research sites, and reduced tug 
charges . 

In the past twenty-five years, ports have become increas- 
ingly crowded, labor and fuel costs have greatly increased 
tug boat charges, and ships have been built that are much 
larger and less maneuverable. Thus lateral thrust units 
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have not only become economical, but they are standard 
equipment for many new-construction ships. 

As more bow thrusters came into use, it became apparent 
that although they work well at zero or very low forward 
speed, their effectiveness drops considerably as the ship's 
speed is increased [2]. In view of the large number of 
these units currently installed, there is suprisingly little 
information available concerning thruster/hull interactions. 

Soon after the first unit was installed in 1955, it was 
learned that the operation of a tunnel-type thrust unit 
causes a low pressure area to form on the ship's hull [2] 
(see Fig. 1) . To find the net resultant thrust, this suc- 
tion force must be integrated over the effected area and 
subtracted from the thrust of the unit. When a thrust unit 
is operated, water surrounding the jet is entrained, as in 
an eductor, causing a low pressure area to be formed on the 
ship's hull in a small area about the jet. In addition, if 
a ship is moving forward, a plume is developed aft of the 
thrust unit and a system of vortices is formed. These vor- 
tices induce a velocity in the water at the ship's hull, 
causing a pressure drop. Additionally, the jet itself acts 
as an obstruction to the flow of water past the ship, caus- 
ing the flow to change speed as it passes, thus causing 
another pressure variation. These induced velocity and 
blockage effects are so severe that a thrust unit's 
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effectiveness can drop as much as fifty percent with a speed 
increase of only three knots [1] . 

The effects of bow thruster jet/hull interactions upon 
a ship motion can be better understood by considering the 
equations of motion for a ship moving in a horizontal plane. 
The general forms of the equations are [3] : 



The linearized forms of these equations are sufficient 
for this analysis and, making the ship's center of gravity 
the coordinate system reference point, they become 
respectively : 



m(u-V(ji) = X Surge equation 

m(v+ui|») = Y Sway equation 



I ib = N 
Z r 



Yaw equation 



-X u (u-u Q ) + (m-x^)u = X 



( 1 ) 



-Y v v + (m-Y^)v - (Y r -mu Q )r - Yd = Y 



( 2 ) 



-N v v - Nd - N r r + (I z -N f )r = N 



(3) 



where: 



m = mass of the ship 



u = velocity in the x-direction (initial value u Q ) 
v = velocity in the y-direction 



t = yaw angle 



r = angular velocity = ip 



I = mass moment of inertia 
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X = excitation forces in the x-direction 

Y = excitation forces in the y-direction 

N = excitation torque 

Dots indicate differentiation with respect to time. 

The subscripts u, u, v, v, r, r, denote differenti- 
ation with respect to the subscript variable. 



The left sides of these equations contain the velocity 
and acceleration-dependent terms such as added mass, added 
inertia, and hydrodynamic drag. The right sides implicitly 
represent the total forces in the x and y directions and the 
total moment about a vertical axis through the center of 
gravity. Although these terms are used to represent various 
excitation forces acting on a ship due to rudders , stabi- 
lizers, propellors, etc., this analysis will deal solely 
with the forces and moments caused by a through-tunnel 
bow-mounted thrust unit. 

In the y-direction> the force due to a bow thruster may 
be written as: 



Y = Y 



t 



+ Y 



s 



where : 



Y 



t 



-m . u . 



jUj, the ideal jet thrust 



Y g = the induced "suction" force on the hull 



m. = the mass flow rate of the jet 



Uj = the jet velocity (positive to starboard) . 
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Inspection of Fig. 2 reveals that -Y g /Y is nearly linear 



ratio range as an example, an approximation of the right 
side of equations (2) and (3) can be made. This portion of 
the curves is of interest because the largest decrease in 
thruster performance takes place here. Since the slope and 
intercept of this line are dependent on hull shape, an 
average performance curve has been plotted in Fig. 3 to 
represent a "typical" hull. A straight line having a slope 
of 5.3 and an intercept of -0.65 is used to approximate the 
curve in the lu^/u^l ratio (hereafter R) of .1 to .3. 

Thus : 

-Y /Y = 5.3R - .65 

u L 

Then since 



By substituting this equation into equation (2) , the 
motion in the y-direction is explicity coupled to the motion 
in the x-direction through the dependency upon R. 

The total turning moment N is the moment produced by 
the force Y acting at a distance x c from the center of gravity. 




Y 



t 




Y_ = m.u.(5.3R -.65) 
s 3 D 



(4) 



and 



Y = -m.u. (1.65 -5. 3R) 
3 3 



.1 * R £ .3 



(5) 
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This characteristic moment arm varies non-linearly as the 
ship's speed changes and is given as: 



x 



Y 



c 



where 



x = the distance from the center of gravity 
to the thrust unit 

x g = the distance from the center of gravity 
to the center of action of the hull 
suction force. 



Chislett and Bjorheden [2] found that x g varies linearly 
with ship's speed. The curve in Fig. 4, which is adapted 
from their data, has a slope of -.7 and an intercept of 1.0. 
Therefore : 



When equations (4) and (5) are substituted for Y g , Y and Y, 
the location of the center of action of the combined force 
system is given by: 




7R) 



R < 1.2 



( 6 ) 



x 



c 



(3.7 1R 2 - 5.76R + 1.65) 

L x 

(1.65 - 5.3R) 



t 



.1 < R < .3 (7) 



and 



N = Yx 



c 
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Thus the yaw equation is coupled to the sway and surge 
equations. Note that equation (7) is only valid for velocity 
ratios between .1 and .3 for this particular fictitious hull. 

Although this analysis has been greatly simplified, it 
demonstrates the complication of the relationships between 
turning moments and ship's speed when a lateral thrust unit 
is used. It can be seen from Fig. 2 that model tests of 
each ship type are necessary in order to accurately predict 
the performance characteristics of a lateral thrust unit. 

With a thorough understanding of the fluid mechanics of bow 
thruster systems, thrust units could be designed to reduce 
the hull suction effect and thus make them more usable under 
a wider range of conditions. An accurate analytical method 
for predicting jet/hull interactions would be useful in 
supplying this information without lengthy and costly model 
tests. Unfortunately, an accurate model that can be super- 
posed with other effects is not available. It is the purpose 
of this thesis to propose an analytical method for predicting 
the jet/hull interactions of a through-tunnel type lateral 
thrust unit. 
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II . BACKGROUND 



As demonstrated in the preceeding section, the ship's 
hull/thruster interactions, when a ship has way on, can 
seriously degrade a bow-thruster's effectiveness. In order 
for naval architects to design a bow thruster installation 
so that it can be used effectively when the ship is moving, 
a complete knowledge of these effects is required. Unfortu- 
nately, only a few waterborne tests have been performed. 

These experiments (Fig. 2) , while clearly demonstrating the 
overall effects on turning ability, do not provide informa- 
tion on the actual pressure distribution on the hull. For 

an incompressible flow, the jet path and the induced flow 

2 2 

are dependent mainly on the momentum flux ratio PjU^ /p^U^ 

and are relatively independent of the Reynolds number [4]. 
Thus an air jet injected into a crossing air flow will induce 
a reaction similar to that of the water jet of a thrust unit 
at the same velocity ratio. There have been a number of 
studies of air jets emitted from a flat plate into a crossing 
air flow in connection with the design of VSTOL (Vertical and 
Short Take Off and Landing) aircraft. Notable studies have 
been made by Bradbury and Woods [4] and Fearn and Weston [5]. 

While the results of these experimental studies are use- 
ful to the naval architect when making an initial design 
intended to improve thruster efficiency, they provide no 
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method for evaluating the resultant configuration. Lengthy 
and expensive model or shipboard tests are required for this 
purpose. An accurate and versatile analytical model of the 
flow interaction effects would provide a quick and inexpen- 
sive means of evaluating new designs, thereby reducing the 
required number of model tests. Some empirical and semi- 
empirical models have been developed. While they can predict 
the pressure distributions on a plate, they cannot be easily 
superposed with other effects such as a secondary jet. 

Because a fluid jet injected into a crossing stream is a 
turbulent, three dimensional, highly non-linear problem, 
empiricism to some degree has been used in all models to 
date. In this study an effort has been made to eliminate or 
at least reduce the empirical inputs to existing analytical 
models while improving their accuracy. 

A review of available models revealed that Wooler's formu- 
lation [6] minimized empiricism, requiring only an experi- 
mentally determined jet trajectory. Hence his model, as 
modified by Waterman [7] , was selected as a basis for further 
investigation . 

In order to visualize the interaction field and determine 
the accuracy of an improved model, once formulated, accurate 
experimental results were required. A review of available 
studies revealed that the pressure contours* for a given 

* 

In this thesis, values of local pressures, P, are expres- 
sed in the form of pressure coefficients, C , where: 

P-P p 

C = — “z 
P J 2P u 

^ CO 00 
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velocity ratio vary considerably from one report to the next, 
as illustrated in Fig. 5 and 6. This difference is especially 
noticeable in the downstream region for values of lUj/U^I 
(hereafter VR) of about 4 , and is present to some degree at 
all velocity ratios. Figure 5 seems to imply that the induced 
flow may not be solely a function of velocity ratio since the 
different studies were not made at the same freestream 
velocities. However, considering studies made with the cross- 
flow velocity nearly equal, as in Fig. 6, it is seen that the 
discrepancies are more probably due to variations in the 
equipment and test procedures used. This serves to point out 
the uncertainty of the available experimental results. The 
results of Fearn and Weston [5] agree well with those of 
other researchers over a broad range of velocity ratios. 
Additionally their report was the most complete of those 
reviewed. It contains not only the usual plots of pressure 
contours, but also raw data, time variation of pressure 
information, and a comparison with other reports. For these 
reasons, the results reported in [5] were selected as the 
basis for determining the accuracy of the analytical model. 

Before attempting to develop an improved model, a thor- 
ough review of Wooler's and Waterman's work was made in 
order to understand their approaches and verify their results. 
A rederivation of these basic formulations will not be pre- 
sented here as it can be found in references [6] and [7] . It 
was noted that Waterman's adoption of a conforming vortex 
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model was a considerable improvement to Wooler's formulation. 
In the conforming vortex model the vortices are constrained 
to lie within the jet rather than being tangential to the jet 
axis as in Wooler's model. By making this alteration, the 
horseshoe vortices are made to conform to the curvature of 
the jet described by the empirical trajectory equation (see 
Fig. 7) : 



5 = B [cosh ( Z/Bd) - 1] (8) 

d 

Where : 

x = the distance down-stream from the jet center, 
d = the jet diameter. 

z = height of the jet axis above the plate. 

2 

B = trajectory parameter = .19(VR) . 

The interference velocity at an arbitrary point on the 
plate is determined by applying the Biot-Savart Law. The 
trailing vortices are approximated as straight line segments. 
As these segments become small, the linear approximation can 
be made to conform to the curvature of the jet. The total 
interference component at a particular point due to the 
entire jet results from the summation of all the elemental 
contributions . 

A comparison of the results of both formulations reveals 
that Wooler's work more closely agrees with the experimental 
findings of Bradbury and Woods [4] than does Waterman's, 
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even though the trailing vortices in Wooler's model do not 
follow the jet trajectory as in the actual case. This result 
is, in part, due to the fact that in Wooler's formulation the 
perturbation velocity at a particular point on the plate is 
caused by vortices acting at different distances from the 
plate. In Waterman's model the interference velocities 
result from the trailing horseshoe vortex contributions act- 
ing at a single radius. Thus the interference velocities 
predicted by Waterman's conforming vortex model are generally 
higher and hence produce lower coefficients of pressure than 
those of Wooler. In addition to modifying Wooler's model, 
and noting studies by Keffer and Baines [8] and Schmidt [9], 
Waterman introduced the concept of an effective origin of the 
vortex sheet at a distance, Z^, above the plate given by 
the furmula: 



Z 

e 



15 



(9) 



Applying this expression to the conforming vortex model. 
Waterman produced the same, if not better far-field agreement 
(Fig. 8) . However, the model lacked accuracy in the near- 
field. 

Waterman further noted that the discrepancies observed in 
the upstream region can be at least partially attributed to 
the blockage of the freestream flow by the jet. It was noted 
that the downstream area is dominated by separation, vortex 
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shedding, and turbulent entrainment. Waterman interpreted 
the blockage effect as that of a circular cylinder in an 
irrotational flow. Using potential flow theory to predict 
the velocity perturbations in the uniform flow caused by a 
doublet centered at the jet origin, Waterman, after some 
numerical experimentation, arrived at an estimate for an 
effective jet diameter given by .96 V vR . The results of 
this model (Fig. 9) are significantly better than Wooler's 
model in the upstream region and far field. This improve- 
ment was accomplished at the expense of introducing two 
additional empirical relationships, the position of the 
effective source and the effective jet diameter. 
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III. IMPROVEMENT IN THE EXISTING MODEL 



Evaluation of Waterman's "fully-corrected conforming vor- 
tex model" (Figs. 9 and 10) revealed that the correlation of 
the analytical model with the experimental data in the area 
just forward of the jet is adequate at each of the velocity 
ratios reported. However, his model is accurate in the far- 
field only for the velocity ratio of 8.0 and in the near- 
field to the sides of the jet for a velocity ratio of 4.0. 

The correlation is poor in the area aft of the jet for all 
velocity ratios. It was apparent that the near and mid- 
field regions, especially aft of the jet, required the 
greatest attention. 

In order to gain a better understanding of the effects of 
the empirical inputs to Waterman's model, it was run a num- 
ber of times using various values of the constants in the 
trajectory, effective source, and effective diameter 
equations. It was learned that small changes in the trajec- 
tory parameter, B in equation (8) , cause little change in 
the overall pressure distribution. The correlation with 
experimental data was improved at some velocity ratios while 

it diverged at others for a particular value of B. A value 
2 

of B = .19(VR) , as before, resulted in the best overall 

correlation of the pressure distribution on the plate. 
Changing the effective diameter of the jet causes 
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major changes in the near-field pressure distribution, but 
has little effect in the far-field. This was expected con- 
sidering the short range of the disturbance caused by a 
doublet in a uniform flow using inviscid relationships. The 
model proved to be particularly sensitive in both the near 
and far-field regions to changes in the position of the 
effective source. Since the quality of the model depends 
so greatly on the effective source correction, the equation 
for it was investigated more closely. 

Equation (9) , the empirical formula for the effective 
source height, was formulated by Schmidt [9] who refers to 
Jordinson [10] and Keffer and Baines [8] for the experimental 
data used as a basis for his work. The Jordinson paper could 
not be obtained for review. The data presented by Keffer 
and Baines, however, does not fit the equation presented by 
Schmidt. The equation, Z e /d = . 8 V VR fits the Keffer and 
Baines data much better as seen in Table 1. 

TABLE I 

COMPARISON OF VALUES OF Z /d 

e / 



VR 


Schmidt 


Experimental 
Keffer & 
Baines 


. 8 > r W 


4 


.53 


1.59 


1.60 


6 


.97 


1.95 


1.93 


8 


1.51 


2.30 


2.26 
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However, Z^/d can vary from .54 V VR to .8 Vvr and still be 

in the range of the heights given in Table I. Using an 

effective source location of Z e /d = . 8 V VR in the Waterman 

model significantly improves the correlation. But, if 

Z /d = . 6 V VR is used, the correlation is much better even 
e 

without the effective diameter correction. Thus, by estab- 
lishing the effective source at . 6 V VR diameters above the 
plate, the full jet diameter can be used in the model and 
the diameter correction is eliminated. See Fig. 11 for an 
example of the improved correlation resulting from this 
change . 

Concentrating on the area aft of the jet, it can be seen 

by comparing Fig. 8 with Fig. 10 that this area is primarily 

affected by the blockage near the plate due to the jet. In 

making multiple runs of Waterman's model and varying the 

parameters of the vortices within the plume, it was confirmed 

that they have little effect on the near-field pressures. 

The far-field to either side of the jet is dominated by the 

vortices within the plume. Therefore, to improve the model's 

accuracy in the near-field, the analysis centered on the jet 

blockage effect. Wu and Wright [11] proposed that a 

fictitious afterbody exists downstream of the jet. This 

afterbody encompasses the highly turbulent, three-dimensional 

flow caused by vortex shedding and entrainment of fluid 

within and downstream of the jet. Outside of this area, the 

fluid streamlines are fairly smooth and close to those of an 

inviscid flow. With this concept in mind, consideration was 
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given to using an elongated shape such as an ellipse or a 
symmetrical foil to approximate the blockage effect. 

In order to establish a probable shape for the afterbody, 
the coefficients of pressure along the corssflow axis from 
the forward edge of the jet to seven diameters aft of the 
jet were plotted using the data of [5] (Fig. 12) . These 
pressure contours resemble those in reference [12] and [13] 
for a symmetrical foil in a uniform flow. Thus a symmetrical 
foil was selected to approximate the jet and the afterbody 
near the plate. 

To determine a length for the foil, it was assumed that 
the foil should end where the coefficient of pressure reaches 
an undisturbed (zero) value. It was not expected that the 
coefficient of pressure should necessarily rise above zero 
aft of the jet due to severe dissipation in the near-wake 
region and the interference effects of the vortices within 
the plume. Indeed the coefficient of pressure leveled-off 
at approximately -0.1 for all available data. Therefore, 

-0.1 was selected as a first approximation of the value of 
the coefficient of pressure used to determine the length of 
the foil. The corresponding distances from the jet center 
to the trailing edge were determined to be 6,2, 3.0, and 2.5 
jet diameters for the velocity ratios of 3.9, 8.0, and 10.0, 
respectively. 

For simplicity, a standard Joukowsky foil was selected 
since it can be easily approximated through mapping methods. 
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Referring to Fig. 13, the Joukowsky transformation, 

Z = 5 + 4 _ 2 

is used to transform the foil in the z-plane to a cylinder 
in the c~plane so that the foil can be treated as a doublet 
in a uniform stream. The perturbation velocities due to the 
presence of the foil are then: 



u 



U 



oo 




U' + a) 

U + D 



- 1 



and 



v 



U 



- J 



h \ u' 


' + a)l 


h') u 


+ 4) J 



Where : 

u = the surface perturbation velocity in the direction 
of the cross-flow due to the presence of the jet 
and afterbody as modeled by the foil. 

v = the surface perturbation velocity perpendicular to 
the cross-flow due to the presence of the jet and 
afterbody as modeled by the foil. 

4 = 1/4 the distance from the forward focus of the 
foil (corresponding to the point (-4,0) in the 
C-plane) to the trailing edge. 

a = the radius of the transformed circle. This term 
is a function of r, the distance from the forward 
focus of the foil to the leading edge. 



26 



Appendix A contains a complete derivation of the above 
transformation . 

The thickness of the foil is a function of both r and 4. 
As described above, the foil length is a function of the 
velocity ratio so that the length l is given by: 



4 = _6_ 
d VR 



(10) 



Initially it was assumed that the forward focus of the 
foil should be at the center of the jet. This makes r equal 
to one-half the jet diameter. When substituted for the 
doublet in the fully corrected model, this method provided 
significant improvement in the downstream area. The up- 
stream correlation was poor, however, since the foil was too 
thick and too blunt at the leading edge. It was necessary 
to adjust the value of r to provide a thinner and more 
streamlined body. Next a foil with its thickness equal to 

one jet diameter was tested. Using the formula, maximum 

2 

thickness/d = 5.2 (m/£ - .7 (m/2.) ) [13], where m is a function 
of r. It was found that if r/d = . 12d/2., where 4 is given by 
equation (10) , the thickness is nearly equal to the jet 
diameter for the entire range of velocity ratios. The foil 
is then moved forward so that its leading edge coincides 
with that of the jet. This configuration (App. B) produced 
the results shown in Figs. 14, 15, and 16 in comparison with 
experimental values [5] for velocity ratios of 4.0, 8.0, and 
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10.0 respectively. Radial and circumferential pressure 
distributions for velocity ratios of 3.9 and 8.0 are plotted 
against the experimental findings of Fearn and Weston [5] 
in Figs. C-l through C-52. 

Attempts were made to use other fictitious shapes such 
as ellipses and ovals as well as distributed source-sink 
combinations. The foil shape, however, produced the best 
results. So the simple relationships for the near-field 
blockage effect given above were retained in the final 
model . 

Reviewing the results, it is seen that for a velocity 
ratio of 8.0, the correlation between this model and the 
experimental data is very close except in the area of 
highly turbulent flow immediately aft (within two diameters) 
of the jet. This is a region of uncertainty for all velocity 
ratios due to the three-dimensional flow caused by entrain- 
ment, vortex shedding, and other highly-viscous wake 
phenomena. Since this area is small, it contributes little 
to the overall thruster/hull interaction. For a velocity 
ratio of 3.9, the correlation is excellent except for the 
mid- and far-field regions aft of the jet centerline. 

This was expected since it was determined earlier that 
the blockage effect of the jet has little effect in the 
far-field. It was observed that if r is such that it causes 
the foil to have a circular shape at the leading edge the 
correlation in the area forward of the jet is improved. 
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However, the agreement deteriorated in the areas to the 
sides and aft of the jet. By varying r/d as a function of 
the location of the point at which the velocity perturbations 
are being calculated, the shape of the foil can be altered 
so that it is nearly circular at the leading edge and becomes 
thin aft of the jet centerline. This configuration worked 
better, but it failed to significantly improve the overall 
pressure distribution resulting from r/d = . 12d/£. This 
observation did lead to the idea of using a half-cylinder 
shape for the jet blockage. That is, the fluid forward of 
jet "sees" a cylinder but once past the jet centerline 
the influence of the cylinder is rapidly diminished. The 
results of this configuration are presented in Fig. 17 from 
which it can be seen that while this is a significant 
improvement over Waterman's model, it is not as accurate 
as the full foil-shape model. It is concluded that further 
improvement of the model must involve a more sophisticated 
treatment of the vortices within the plume as well as the 
blockage effect of the jet. 
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IV. CONCLUSIONS 



As evidenced by Figs. 14, 15, 16, and C-l through C-52 
this analysis adequately predicts the pressure distribution 
on a flat plate from which a jet is emitted into a crossflow. 
Application of the model is limited, however, to velocity 
ratios VR = 3.9 or greater. This is due to the method by 
which the length of the fictitious afterbody is determined. 
Above velocity ratios of about 4, this length decreases for 
increasing VR, but plots of the coefficients of pressure 
along the crossflow axis similar to Fig. 12 indicate that 
the afterbody length should decrease also for values of VR 
less than about 4. This is consistent with the observations 
of Fearn and Weston. They reported that an anomaly in the 
pressure distribution in the wake occurs at a velocity ratio 
of 3.9. The implication is that the rate of pressure recovery 
in the wake region is a minimum for a velocity ratio of about 
3.9 [5]. 

No attempt has been made to accurately predict the pres- 
sure distribution in the highly turbulent wake of the jet. 

This area, having a width of about one jet diameter, extends 
downstream from the jet along the crossflow axis. The pre- 
dicted values of the plate pressures reported here cannot be 
relied upon since most of this area is contained within the 
fictitious afterbody. 
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Figure 18 shows the results of this analysis plotted as 
the ratio of the induced side force to the jet thrust, 

-Y /Y^, versus the velocity ratio VR. These results are 
compared with NPL pressure measurements on a flat plate and 
the model tests of English [1] and Chislett and Bjorheden [2] . 
The slope of the predicted -Y /Y curve is approximately 5.5, 
whereas that of the experimental data is about 5.3 for both 
the flat plate and the model tests. While this analysis 
predicts somewhat larger values for the overall (integrated) 
interference effects, the small difference between these 
slopes indicates that the nature of the relationship between 
the interference effects and the velocity ratio predicted by 
the analytical model is correct. Fearn and Weston did not 
report overall interference effects so the correlation between 
their work and the NPL experiments is not known. Since the 
assumptions concerning the shape and size of the afterbody 
were based on Fearn and Weston's data, the degree to which 
this analysis is in error cannot be determined from Fig. 18 
alone. Even with these deficiencies, it is felt that this 
analytical model adequately predicts the pressure distribu- 
tion on a flat plate and can be used as a first approximation 
for design purposes. 
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V. RECOMMENDATIONS FOR FURTHER STUDY 



In Chapter III it was noted that afterbody shapes other 
than a symmetric foil have been tried without significant 
improvement. This does not mean that this approach showed 
no promise for possible improvement in the model, but only 
that the actual shapes tested were not the best. It is felt 
that this model can be improved by changing the shape used 
for the jet and afterbody. 

As reported earlier, the crossflow fluid particles near 
the plate apparently see a cylindrical shape when approaching 
the jet. As they pass the jet, some particles are entrained, 
some join the turbulent flow within the wake, and, further 
out, they pass as if the jet and wake together are a solid 
object. Just downstream of the jet, along the crossflow 
axis, there is an area of high entrainment. This causes a 
high local pressure drop out to about one jet diameter for 
VR = 8. Pressure recovery is rapid past this point. This 
leads to the idea of using a shape such as: 




This shape has a cylindrical leading edge, the maximum 
thickness is just downstream of the jet center, and then it 
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narrows quickly to a relatively thin tail. The total length 
is the same as that for the symmetric foil. 

It is expected that this shape can be formed using the 
Joukowsky transform of the existing model by varying r as a 
function of the angular position of the particular point at 
which the perturbation velocities are being calculated. 
Although a functional relationship which will work over a 
range of velocity ratios has not yet been found, it is pos- 
sible that one exists. This idea might also be pursued 
using panel (distributed singularity) methods. 

Another area of uncertainty during the course of this 

study was the empirical jet trajectory equation. Reports of 

experiments to determine the jet axis provide a wide range 

of possible trajectories. If the trajectory equation used 

in this model is accepted as correct, the trajectory con- 

2 

stant B = . 19VR still may not be accurate over a wide range 
of velocity ratios. In particular, it is suspected that it 
may be in error for velocity ratios near 4 where an anomaly 
in the plate pressure distributions occurs. The cause of 
this anomaly is not clear, but it is suspected that this is 
the point where the jet momentum is so low compared to that 
of the crossflow that the jet reattaches to the plate within 
a few jet diameters. If so, the jet trajectory would be 
quite different than that which is currently modeled. 

Finally, the entrainment of the crossflow fluid into the 
jet is implicity accounted -for in the formulation of the 
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horseshoe vortices. However, the pressure drop associated 
with the jet entraining the surrounding fluid is not explic- 
itly determined. This model predicts a zero pressure distri- 
bution at very high velocity ratios (zero ship speed) . 

Further improvement of this model in accounting for entrain- 
ment is needed if it is to be used to predict the performance 
of bow- thrusters in ships with no way on. 
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Hull low pressure areas at various speeds. 




Figure 2. Hull suction effects vs. R. 
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Figure 3. -Y /Y fc vs - R for a "typical" hull. 
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Figure 4. Center of action of the suction force. [ 
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Figure 5. Comparison of experimental findings, VR= 
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Fiyure 6. Comparison of experimental findings. 



Figure 7. Jet axis and conforming vortex loop system. [7] 
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Figure 9. Comparison of pressure coefficient contours, experiment 
vs. fully corrected conforming vortex model, VR=8. 
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Figure 10. Comparison of pressure coefficient contours, experiment 
vs. fully corrected conforming vortex model, VR=4. 
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Figure 12. Plots of coefficients of pressure along the free- 
stream axis used to determine the length of the 
foil. Data from [5]. 
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Figure 13. The Joukowsky transformation. 
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Present model 
Fearn & Weston [5] 
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Figure 15. Comparison of pressure coefficient contours, experiment 
vs. present analytical model, VR=8. 
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Figure 16. Comparison of pressure coefficient contours, experiment 
vs. present analytical model, VR=10. 
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Figure 17. Conforming vortex model with a half cylinder for 
blockage correction vs. experiment, VR=8. 



NPL 

flat plate 




Figure 18. Comparison of theoretical values of -Y g /Y t with 
experimental findings. [1] 
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APPENDIX A 



DERIVATION OF THE JOUKOWSKY TRANSFORMATION 



2 

The Joukowsky transformation, Z = c + l /?, is used to 



transform a foil in the z-plane to a cylinder in the c-plane 
so that the foil can be treated as a doublet in a uniform 
flow. The terms used in this transformation are: 

Z - the complex position vector of a point in the 
Z-plane Z = x + iy 

C - The complex position vector of a point in the 
C-plane. c = £ + in 

l - A positive real constant equal to 1/4 the distance 
between the foci of the foil. 

a - The diameter of the cylinder in the s' -plane. 

m - The distance between the n and n ' axes. This is a 
function of r. 

r - The distance from the forward focus to the leading 
edge of the foil. That is r = |c - 2l\ where -c is 
the abscissa of the leading edge. 

C ' -plane - The linearly transformed t;-plane used to give 



From the Joukowsky transformation and Fig. 13 it is seen 
that : 



the foil thickness, c' = e,' + in' 




(Al) 
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m 



.5 (A 



l) 



(A2 ) 



C ' = C - (~m) = £ + m 

The complex potential of a cylinder in a unifor 
is given by : 



tO 



(c') = U ao (c' + a/?') 



The complex velocity is then found as follows 



Uoo + u - iv = 



dto 


dto 


dc’ 


dZ 


dc' 


d£ 


dto 


dto 


dc ' 


dZ 


dc' 


dc 



dc 

dZ 



dZ 

dc 



then substituting 



dto 

dZ 



-°-(W 



( £ ' 2 - a 2 ) 
(£ 2 - * 2 ) 



and since £ ' - a = £ - 2. 



dto 

dZ 



= U 






( g * + a) 
(£ + 4 ) 



(A3) 

stream 



(A4) 



(A5 ) 



(A6) 
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G|< G|C 



Then 




U' + a) 1 

(?' + 4) J 



1 



(A7) 




(C + a)l 

U + J 



(A8) 
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APPENDIX B 



COMPUTER PROGRAM 



C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r 

C 

c 

c 

c 

c 

r 



XX XX XX XXX XXX xxxxxxxxxxxxxxxx 
X X 

X THIS PROGRAM SIMULATES A X 
X FLUID JET INJECT=0 IMTn x 
X A UNIFORM CROSSFLOW AND X 
X FURTHER OE t FRMINES TH c X 
X PRESSURE 0 ISTRIBUTION ON X 
X THE SURFACE r-ROM WHICH X 
X THE JET IS EJECTED X 

X X 

XX XX XX XXXXXX XXXXXXXXXXXXXXXX 



CEFINITICN OF TERMS 

NR AD - NUMBER OF RADIAL POSITIONS DM SURFACE 

( I NPUT) 

NANG - NUMBER OF ANG'JIAR POSITIONS PM SURFACE 

( I NPUT ) 



NP 0 I NT - NUMBER OF ARBITRARY POINTS ON S U R c -A C E 
(MR AQ*NA NG) 

UJ - JET VELOCITY ( TNPUT) 

UM - MAIN STREAM VELOCITY (INPUT) 

0 - ACTUAL OIA , -’FT = R OF JET nDjcifc (THRU!) 



- NON— D TM C NSI DMA I RADIAL D OS! T TPN OF SUE F A C c 
POINTS (DIMENSIONED AT LEAST NR AD ) (INPUT) 



TH FT A 



XP 

YP 

X 
Z 

Z I NC 

OELTAZ 

HKSAV 

XI 

ZETA 

ETA 



AL °H A 



ANGULAR P 0 SITIP\| nc SUPFAC r ROTM^S 
(OIMFMSIOVFD AT LEAST NAMQ) 

X-°OS T T I ON OF SURFACE POINTS 
(DIMENSIONED AT LEAST NPOINT) 

Y -°OS I T I ON OF SURFACE POINTS 
( DI M C NS TO MFO AT LEAST NPOINT) 

HORIZONTAL COORDINATE OF J=T AXIS 

VERTICAL COORDINATE 0 F JET AXIS 

NON-OI‘*ENS ! CNAL INCREMENTAL STS° SIZE for z 

z-position OF E c F ECTI VE/VIR TUAL SOURCE 

CUMULATIVE VCRTICITY ALONG THE J C T 

NATUFAL COORDINATE TANGENT TO JET 

NATURAL COORDINATE NORMAL TO J C T 

NATURAL COO POT NAT c PERPENDI CULAP T 0 XI 
AND ZETA 

ANGLE BETWFFN X-AXIS AMD XI-AXIS 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



PA P E A - INCREMENTAL PLATE AREA 
TDTFCR - TOTAL RESULTANT FORCE ON PLAT 



X8 AR 



X-RQS I T I ON OF RESULTANT FO^CP 



YT 



NORMAL I ZED INDUCED FORCE 



FLL 



1/4 THE DISTANCE BETWEEN TH c FOCI 
OF THE FOIL USED FOR THE BLOCKAGE EFFECT 



ra nj 



ZC? 



THE DISTANCE FROM THE FORWARD FOCUS OF THE 
FOIL TO THF LEADING EDGE 

THE COMPLEX VE LOC I T Y INDUCED IN TH = UNIFORM 
FLOW DUE TO TH E PRESENCE OF THF FOIL 



IMPLICIT REAL *9 (A-H,0-Z) 

01 MENS I ON U( 50 ,50) , V ( 50 , SO ) , W ( 50 , 5 0 ) , R ( SO ) , THET A { 50 ) , 
X VFL (450) . C P ( 5 0 . 5 0 ) 

01 MF NS I CM X P ( 4 5 0 ) ,Y° (450), Z°(450) 

COMPLEX ZED* 16 ,ZST*16»ZC°* 16 ,ZSTP*15 , ZR*1 6 , 
*ZEL*L6,Z4*16,ZHA*16, ZM*16 
C DATA INPUT 

RE AD ( 5 , 60 ) NP AO 
READ (5 , 60 ) NANG 
6 C F0RMATU4) 

NP 01 NT = N R A D* N A NG 
C ZERO OUT STORAGE MATRICES 
DO 3 0 I = l, MR AD 
DO 2 0 J = l, NANG 
U( I , J) =0.0 
V< I , J) = 0.0 
C° ( I ♦ J ) = 0.0 
20 CDNTTNU c 

30 CONTINUE 

RE AO (5, 40 ) (JJ.UMtO 
WR I T E ( 6 ,5 0 ) UJ ,UM , 0 
C VELOCITY RAT to 
VR =(UJ/')M) 

READ (5, 70) ZINC 

R= AD ( 5 , 70 ) (P ( I ) , 1 = 1, NRAD ) 

70 FORMAT (FLO .5 ) 

C SET CONSTANTS US c l> IN FCIL BLOCKAGE CALCULATIONS 
31 =1.0 
8HA = 0.5 
84 = 4.0 

C SET p 0 I L PARAMETERS 
ELL = o.O/VR 
R A 0 J = . 1 2 /E LL 

ALAM = .5*(2.*ELH-R4DJ+0SQRT( ( 2.*ELL+RA0J )**2- 
*4. 0*ELL**2 ) ) 

EM = 0. 5* ( AL AM - ELL ) 

C POSITION PF EFFECTIVE SOU PC E 
D c LT az = . 6*DSQ!*T( V R ) 

WR I T E ( 6 , 3 0 ) DFLTAZ 
WRITE (6, 90) ZINC 
40 FORMAT ( 3F 3.4 ) 

50 FORMAT( *0* ,2X, ' JET VELOCITY = ', F 8 . 4 , 2 X , ' UN I p 0 D 1 , 
X'VELCCITY = ' ,F3.4,2X, ' JET n I A M E T E R = »,F3.4) 

80 FORMAT! 'O' ,2X, ‘EFFECTIVE SOURCE IS • , F5 . 2 , 2X , 
X*dTAVFTERS A BOV - TH C PL * T E ' ) 

90 FORMAT {‘0',2X.’Z STEP SIZE = ' , c 5 . 3 ) 

C CFFINITION of Pt 

PI = 2.0* A RS IN (1 .000 00000) 

C TRAJECTORY CONSTANT 
8 = 0.19*VR**P 
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c 



c 

c 

c 



c 

c 

c 

c 

c 

c 

c 

c 



CALCULATIC.N OF INTERFE C PNC 5 VELOCITY 
L = 0 

CALC ON = OFLOAT ( MANS ) - 1.0 
INCREMENT RADIUS (INPUT VALUES) 

DO I AO I = I, NR AD 
PC AL C = » ( I ) 

INCREMENT THETA (5 DEGREE INCREMENTS) 
DO 130 J = I, NANG 
K = J-I 

THETA(J) = DFLO AT ( K ) * <°t/CALCCN) 
ARAD = ALAN - EM 
SHIFT TO CARTESIAN COORDINATES 
t = L+l 

XP(L) = R CALC* DC OS (THETA( J ) ) 

Y° ( L ) = PCALC*DSIN(THETA( J ) ) 

HKSAV = 0.0 
12 = 0.00001 
100 Zi = 12 

12 = Zl + ZINC 

Zl l = Zl + oeltaz 

122 = Z2 + OELTAZ 

XI = S*(DCOSH( Z 1 / 3 ) - 1.0) 

X2 = B* (OCOSHI Z2/B)- 1.0 ) 

ALFONF = OATANI l.O/D SINH( Z 1/3) ) 
ALFTWO = D AT AN ( 1 . 0 / D S IN H ( Z 2 / 3 ) ) 
ALF3AR = DATAN ( ( Z2-Z 1 ) / I X2 -X I ) ) 

SINAI = OS I N ( At. FON F ) 

SI NA 2 = 0 S I N I A L c T W O ) 

SINAB = DS IN( AL C BAR) 



COSAI = PC OS (A L c ON E ) 

CD S A 2 = ncns ( A L c TWO) 

CO S A 8 = ncns ( ALPPAR) 

NATURAL COORDINATES IN TERMS OF CARTESIAN COORDINATES 
XIONE = -(Zll*S?NAl + ( X l- XP { L ) ) vCOS A 1 ) 

XI TWO = - ( Z22*S I NAP. + (X2-X°(l ))*CRSA2) 

XI BA R 1 = — ( Z 1. 1 * S I N A R -t- ( <1 - XP { L ) ) *CD S A B ) 

XI3AR2 = — ( Z22*S IN A6 + ( X 2 - X P ( L ) ) *CD S A 8 ) 

ETA = YP( L ) 



ZETA1 = (Zll*COSAl 
ZE TA 2 = ( Z22*C0S \2 
Z C TA3 = ( Z1I*C0SAB 
GEGM=TRIC CONSTANTS 



(Xl-X P { L ) ) *S INA l ) 
( X 2 - X P ( L ) ) * S I N A 2 ) 
( X 1 - X P ( L ) ) *S IN A B } 



ET PL = C T A + 0.5 
FT mi = ETA - 0.5 
SA = XI ONE** 2 
SB = X I TW 0**2 
SC = XI C »ARIR:«2 

50 = X I BAR 2**2 

SE = Z-TA1**2 

SE = Z C TA2**2 

SG = Z E T AB**2 

SH = C T Pt **2 

51 = PTMI**2 

DA = DSQPTISA + S c + SH ) 

DB = D SCR T ( S A + SE + SI ) 

DC = DSCPT(SB + SF + SH ) 

DE = DSCRT (SB + SE + SI ) 

DF = DSCRT (SO + SG -t- SI) 
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DG = 


0 SCR T ( SC +• 


SG 4- SI ) 




DM = 


DSCPTISD ♦ 


SG + SH) 




ni = 


OSCRTISC + 


SG + SH) 




CONI 


= (ETRL/DA 


- ETMI/OB ) / ( SA * S c 


) 


COM2 


= (ETPL/DC 


- ET mi/df ) / ( SB 4- SF 


) 


CON 3 


= ( X I B A P 2 / ' 


AF - XIBAR1 /PG) /( SG 


4- 


C3N4 


= (XI3AP2/04 - XI6AR1 /0I)/( SG 


+ 


Ul = 


ZETA 1*C0N l 






U2 = 


Z FT A2*C0N 2 






V3 = 


ZFTA3*CGN3 






V4 = 


Z ETA 6*C ON 4 






W1 = 


X IC.NF*C0N1 






W2 = 


XT T 'JO*C0N2 






W3 = 


F T M j * C. 0 N 3 






W4 = 


E T D L*C.0N4 







C INCREMENTAL VQRTICITY TERM 

FF = ( ( TF X P ( Z2/3)-0EXP( Zl/ 3 ) )/( I .OOEX? ( ( Z2+Z1) /3 ) ) ) 
HK = (L. 0/3.0)* VR**2 * DATAN(FF) 

HK = HK + HKSAV 
HKSAV » HK 

C INCREMENTAL I NT ER F c R E NC E VELOCITY COMPONENTS 
IIP = 0.0 

UP = HK* ( ( — W 1*S l N A 1 + W 2* S I N A 2 + ( -w3+k 4 ) *S I N A8+U l*C 0 S A 1 ) 
X - { U2*CDSA2 ) ) 

VP = 0.0 

VP = HK* ( V3 - V 4 ) 

U( I , J) = UII.J) + HP*2. 0 
V ( I f J ) = V ( I , J ) * VP *2.0 
C TEST FTP CONVERGENCE 

TEST = DA b S{ UP) + OABS( V° ) 

IP ( TEST. LE. 0.000001 ) GO TO 110 
GO to ICO 
110 CONTINUE 

C CALCULATION OP BLOCKAGE VELOCITY T F RM 
C TPANSEORm COO°D INA^PS TO Z PLANE 
8 3 A Z - X P ( L ) — 2 . 0 * E L L + (.5-RADJ) 

IF ( AZ. EC. 0 .0 . AND . Y°< L) . =0. 0.0) GO TO 115 
C CHANGE VARIABLES TO COMPLEX NUMBERS 
3Z = Y?(l.) 

BB = 0 . 0 

Z C Q = OC.MPI.X (A Z » 3Z ) 

ZR = OCMPL X( ARAD, BB) 

ZEL = OCMPLXIELL.BB) 

Z4 = OC VO L X { 34 » 33 ) 

ZM = DC V P LX ( EM , 33 ) 

Z4A = DCMPLX (BHA , BB) 

IF ( AZ .L T .0 .0 )G0 TO 8 9 

C TRANSFORM PQSI t ION VECTOR IN Z PLANE TO Z FT A PLANE 
ZST = ZHA* ( ZED + C OS GRT ( Z5 D**? - Z4*ZE L **2 ) ) 

GO TO 91 

89 ZST = ZHA* ( ZED - CDS CRT ( ZE D**2-Z4*EI.L**2 ) ) 

C TRANS FORM POSITION VE.CTORS FORM ZETA P L A N F TO 
C ZETA PRIME PLANE 
91 ZS T° = ZST + ZM 

C S FT VELOCITIES TO ZERO WHEN WITHIN THE FOIL 
XF. 0 = DR C AL ( ZS T° ) 

YE D = niMAG(ZSTP) 

pp = n SCR T ( X Fn ** 2 + YFn**2 ) 

IF { PR .LF. ( ARAD/1A. 0) )G0 TO 05 
C CALCULATE the DERIVATIVE OF THE COMPLEX POTENTIAL 
95 ZCP= (( ZS T /ZSTP )**2 )* (ZST° * ARAO)/(ZST 4- r -\_L ) 

C EXTRACT THE HPRiz. A^O VEP t . P^RTABATHN V C L°CITIES 
UC = D p EAL ( 7 CP ) - 1.0 
VC = -CIMAGIZCP) 

GO TC 119 
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115 l C = 0 . C 
VC = 0.0 

119 U( T , J) = U( I , J ) * UC 

V( I . J) = V ( I . J ) + VC 

C CALCULATION OF PRESSURE COEFFICIENT 

120 C°(I,J) = -2.0*(U(I,J) + .5*') ( I , J )**2 + . 5*V( I , J )**2) 

C CHANGE DESIGNATION OF VARIABLE 

ZP (l ) = CP ( I ,J J 

VEL(L) = DSQRT((U( I, J) + 1.0 )**2 + V(I,J)**2) 

130 CONTINUF 
140 CONTINUE 

C CALCULATION OF INCREMENTAL FORCE 
TOTFOR = 0.0 
TYMOM = 0.0 
MARE A = NR AD - 1 
00 170 1=1, NAPE A 
. K = I + 1 
PI = R ( I ) 
prj — P ( K ) 

C CALCULATION OF INCREMENTAL AREA 
AR FA = PI* I PC*-? - PI**?) 

PAPEA = A P EA / ( 2 . 0*C A LCON) 

DO 160 J = 1 ,N ANG 

C CALC'JLATICN OF AVERAGE PRESSURE COEFFICIENT 
CPAVG = (CP(ItJ) * CP(K,J ) J/2.0 
C CALCULATION nP FF FEC T IV E FORCE 
FT PCS = ?APEA*CPAVG 
IF(J.FO.l) GO TO 150 
IF ( J.EC.MAf'G) GO TO 150 
FORCE = 2 . 0* FO RO c 
C CALC'JLATICN OF MOMENT ACM 
150 XM = ((RO + R I ) / 2 . 0 ) *0C OR ( TH FT A ( J ) ) 

C CALCULATION OF MO M ENT ABOUT V -AXIS 
VMOM = Fop.CF*XM 
C SUMMATION OF total c ORC E 

TOTFOR = TOTFOP + FORCE 

c SL’ MM ATirN c p Total moment 

TY MOM = TYMOM +• YMOM 

16 0 CO NT INI IF 
170 CONTINUE 

C CALCULATION OF C C NTER CF ACTION 



C 


X-3AR = TYMCM/TQTFOR 
CALCU I.A TI C A OF NORMALIZED INDUCED 


FORCE : Y/t 


c. 


VRINV = 1.0/ Vo 
YT = 2 .0*TCTF0R*VRINV**2/P I 
DATA OUTPUT 




LdO 


WR ITF( 6,230) TOTFOR, XO AP 
WR I T F ( 6 ,240) YT, VRINV 
on ??0 I = l , NANG 




190 


ANGLE = T H C T A ( I )*180 .O/PI 

WRITE (6,190) ANGLE 

c 0 RM AT ( *0 ' , 2X » 'THETA = *,F6.2, 


8X , • R ' , 9 X , 1 C R ' , 9 X , ' U ' , 



X 9 X , • V • ) 

DO 2 10 J = 1 ,NP AD 

WP IT P( 6,2 00) R ( J I ,CP ( J, I ) , U< J, I) ,V( J , I ) 

200 FO F M AT ( • • , 23X , F6.2, 2X, F8 . 2, 3X , F6.2, 4X, f=6.2) 

210 CUNT IN UF 

22 0 CONTINUE 

230 RORMAT( »0« ,2X, 'CESULTANT FORC^ OF • , l X , «=7 . 2 , IX , 
X'ACTS AT X = • , R 6 . 2 ) 

240 FO F. M AT ( • 0 ' , 2X , * Y/T = • , F6 . 2 , IX , ' FOP ijm /U J = «,F6.3) 

WC I T c ( 7 , 2 5 G ) ( XP ( I ) , Y P ( I ) , Z d ( I ) , I = 1 , N PO I NT ) 

250 FORMAT ( 3F20. 4) 

270 ST OD 
END 
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COMPARISON OF PRESENT MODEL WITH 
EXPERIMENTAL RESULTS 
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RADIAL DISTANCE (X/D) 
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RnDIflL DISTRNCE (X/D) 
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RADIAL DISTANCE (X/D) 
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RADIAL DISTANCE (X/D) 
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RRDIRL DISTRNCE (X/D) 
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Figure C-16. C ' s along circumf rential arcs, VR = 8, radius = .58 dia. 



o 



I— 

cn 

LU 



O 

in 

r- z iu 

• CO LJ 







O 

LO 



o 

• 

rr 

i 



CD 

LU 

O 



CL 

I— 

LU 

LL 

1— 



3ynSS3Ud 30 1N3I3I333Q:) 



77 



Figure C-17. C ' s along circumfrential arcs, VR = 8, radius = .75 dia. 
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Figure C-18. C ' s along circumf rential arcs, VR = 8, radius = 1 dia. 
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Figure C-19. C 's along circumf rential arcs, VR = 8 , radius = 1.4 dia. 
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Figure C-20. C 's along circumf rential arcs, VR = 8, radius = 1.6 dia. 
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Figure C-21. C ' s along circumf rential arcs, VR = 8, radius = 2 dia. 
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Figure C-22. C ' s along circumfrential arcs, VR = 8, radius = 2.5 dia. 
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Figure C-23. C 's along circumfrential arcs, VR = 8, radius = 3 dia. 
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Figure C-24. C 's along circumf rential arcs, VR = 8, radius 3.5 dia. 
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Figure C-25. C 1 s along circumf rential arcs, VR = 8, radius = 4 dia. 
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Figure C-26. C 's along circumf rential arcs, VR = 8, radius = 5.5 dia. 
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Figure C-42. C ' s along cricumf rential arcs, VR = 3.9, radius = .58 dia. 
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Figure C-43. C 's along circumf rential arcs, VR = 3.9, radius = .75 dia. 
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Figure C-44. C 's along circumf rential arcs, VR = 3.9, radius = 1 dia. 
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Figure C-45. C 's along circuiafrential arcs, VR = 3.9, radius = 1.4 dia. 
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Figure C-46. C 's along circumf rential arcs, VR = 3.9, radius - 1.6 dia. 
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Figure C-47. C 's along circurof rential arcs, VR = 3.9, radius - 2 dia. 
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Figure C-48. C 's along circumf rential arcs, VR = 3.9, radius = 2.5 dia. 
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Figure C-49. C 's along circumf rential arcs, VR = 3.9, radius = 3 dia. 
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Figure C-50 C 's along circumf rential arcs, VR = 3.9, radius = 3.5 dia. 
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Figure C-51. C 's along circumf rential arcs, VR = 3.9, radius = 4 dia. 
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Figure C-52. C 's along circumf rential arcs, VR = 3.9, radius = 5.5 dia. 
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